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Abstract 



We study J/ip dissociation processes in hadronic environments. The validity 
of a thermal meson gas ansatz is tested by confronting it with an alternative, 
nonequilibrium scenario. Heavy ion collisions are simulated in the frame- 
work of the microscopic transport model UrQMD, taking into account the 
production of charmonium states through hard parton-parton interactions 
and subsequent rescattering with hadrons. The thermal gas and microscopic 
transport scenarios are shown to be very dissimilar. Estimates of J/ifi survival 
probabilities based on thermal models of comover interactions in heavy ion 
collisions are therefore not reliable. 
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I. INTRODUCTION 



The suppression of heavy quarkonia has long been proposed as a unique signature for 
the transient creation of a quark-gluon plasma in heavy ion collisions [|J. Preliminary 
reports on J/ip and ip' production cross sections in Pb+Pb reactions at SPS energies 
in fact seem to indicate an anomalous charmonium suppression compared to pA and lighter 
AB interactions. It has been demonstrated that the framework of Glauber theory with 
some comover interactions fails to describe all the data while assuming the onset of 
deconfmement at a certain interaction density leads to a better agreement with experiment. 
Starting from a similar ansatz, however, it was claimed that hadronic mechanisms alone 
could account for the observed suppression probabilities || . While in these studies universal 
comover absorption cross sections are parameters, there have been attempts to predict the 
actual comover cross sections for certain reaction channels and to estimate the resulting J/ ip 
survival probabilities []7]||. From these studies one learns that the dissociation cross sections 
may depend strongly on the scattering energy as well as the meson species interacting with 
the charmonium state. Unfortunately only thermal absorption rates have been calculated 
on the basis of the predicted cross sections. 

It is the goal of the present work to investigate, on the basis of a hadronic transport model, 
if the assumption of a thermal hadronic environment can be justified. This is analogous to 
studies of the nonequilibrium processes contributing to strangeness production || versus 
thermal model calculations |TIJ. Scattering rates and collision energies in the microscopic 
simulation are compared to the ideal hadron gas scenario. We also address the question of 
whether only the consideration of light mesons (tt and p) suffices to account for the major 
part of all possible comover absorption processes. 



II. THE MODEL 

The production and subsequent absorption of charmonium states in relativistic heavy 
ion collisions are the result of an intricate interplay between hard and soft processes. Since 
the interactions are nonperturbative the calculation of their dynamics from first principles is 
prohibitive. Thus phenomenological models which account for both hadronic and partonic 
aspects are needed. We simulate charmonium production and absorption microscopically. 

The charmonium states are produced exclusively in hard parton-parton interactions for 
which a simple Glauber-type picture of nucleus-nucleus reactions applies. This is supported 
by the fact that the experimentally measured dimuon continuum in the corresponding mass 
range can be well described by perturbative calculations of the Drell-Yan process^. In 
particular, a linear scaling of the Drell-Yan cross section with atomic mass A is observed 
in pA interactions [II]. Initial state interactions are assumed to be negligible. We apply 
perturbative QCD to nucleus-nucleus collisions without nuclear modifications of the parton 



distribution functions. (See Refs. [|T3,[y| f° r the effects of such modifications 



The production of cc pairs which eventually form the quarkonium state is, however, dominated 
by gluon-gluon scatterings. 
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Within perturbative QCD, the hard production process is completely decoupled from 
other possible interactions such as the creation and reinteraction of soft quarks and gluons 
and their rearrangement into color singlet states. In particular, the hard production process 
is not linked to soft hadron production, baryon stopping and hadronic rescattering. 

The dissociation of charmonium states is assumed to be due to hadronic interactions. 
The space-time evolution of mesons and baryons, governed by soft physics, can be treated 
by transport theory on the hadronic level. The kinetics of each interaction is determined 
by appropriate hadronic scattering cross sections. In particular, charmonium dissociation 
can be studied by introducing cross sections for their absorption through interactions with 
baryons and mesons. In this work, we consider dissociation by comoving mesons only[|. 

Technically, we generate a space-time distribution of charmonium production points for 
any impact parameter b by microscopically simulating Glauber-type nucleus-nucleus colli- 
sions assuming that the nuclei, i.e. sets of appropriately initialized nucleons, pass through 
each other on straight line trajectories without energy loss. At the space-time point of each 
nucleon-nucleon encounter, a charmonium state has a fixed production cross section. 

To calculate the rescattering of charmonium states we generate a full hadronic cascade 



simulation with the Ultrarelativistic Quantum Molecular Dynamics, UrQMD, model |15 
which takes stopping and particle production into account. The charmonia are inserted 
into the evolving hadronic environment at the appropriate space-time points according to 
the Glauber simulation of the nucleus-nucleus collision. The momenta of the charmonium 



states are assigned from the parametrization [16 



E ^T71 ~ (! -^) 3 - 55 exp(-p T 2.08GeV- 1 ) . 

In principle, the total energy of the produced cc state has to be subtracted from the overall re- 
action. However, it is difficult to do so consistently since simulations of the charmonium pro- 
duction and rescattering processes are formulated within different conceptual frameworks. 
Therefore we allow for the small violation of energy-momentum conservation imposed by 
this procedure. In Pb+Pb collisions at 160 GeV/nucleon, = Ay/s NN ~ 3600 GeV and 
with a typical charmonium energy of E c5 m 4 GeV, E c5 /y/s ~ 0.1% . 

In the simulations we have used fixed dissociation cross sections, ctj/^m = 2crj° t ^ 7V /3 
according to the additive quark model, for any meson species M. The value of the total 
J/ipN absorption cross section, taken from Ref. |I7|] , is cf Xj jJ , ^ )N = 3.62 mb. Since we do 
not consider elastic interactions, any collision above the kinematic threshold leads to the 
dissociation of the J/ip so that we use the two terms synonymously. For this study, we 
have neglected the energy dependence of the dissociation cross sections which may be quite 
strong close to threshold. Also, we have not considered a possible (eigen)time dependence 
of the actual dissociation cross sections which may evolve from a small initial cross section 
at the production point to its asymptotic value (see 0]). We want to focus on the question 



2 The calculation of nuclear absorption within our model framework shows interesting deviations 
from the standard Glauber-like prescription. Due to nuclear stopping in the UrQMD simulation, 
incident nucleons are slightly retarded which leads to an increase of the effective path length of 
J/ip's in nuclear matter. This will be discussed in 
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of whether or not the J/^-comover reactions in a heavy ion collision can be reasonably 
approximated by a thermal hadron gas interacting with the J/ip. The qualitative answer to 
this question is rather insensitive to the actual parametrization of the cross section. 

Note that we do include the formation time of the comovers (on average, 7> ~ 1 fm/c). 
Particles produced in a string fragmentation process are not allowed to interact with each 
other or with a J/ip within their formation timef^. 

We assume E t h = 2mo = 3730 MeV is the kinematic threshold for J/ip- meson dissoci- 
ation processes. In the thermal model, a fixed p mass of 770 MeV is used, therefore the 
minimum J/ip p collision energy is rrij/^p + m p = 3870 MeV > E t ^. In the UrQMD model 
however, the finite width of the p meson, T = 151 MeV, is taken into account. The p meson 
mass distribution is selected from the Breit-Wigner resonance formula. Therefore, in the mi- 
croscopic model some of the J/ip p interactions are below threshold while no correpsonding 
threshold exists for the thermal model. 

The results in the equilibrium scenario are easily obtained. Consider an ideal gas of 
J/f/Vs and one meson species M (n or p). We calculate the rate R of dissociation processes 
for an arbitrarily chosen J/ip in the rest frame of the J/ip: 

R = 9M J d 3 p M f(p M )j(pM)(Tj/ipM, 

where the spin and isospin degeneracy, gM, is 3 for 7r's and 9 for p's. The mesons are 
distributed according to the Bose-Einstein distribution, 

f(PM) = -^(exp(£ M /T) " I)" 1 



where E M = ym 2 M + p M is the total energy, pm is the momentum of mesons in the rest 
frame of the J/ip and j{pm) = Pm/Em is the meson flux. The dissociation cross section has 
a fixed value, (Tjmm — 2.41 mb. The above integral can be rewritten in terms of the center 
of mass collision energy, E crn , 

R= f dE cm ^p M ^^-(exp(E M /T) - \y x oj^ M = [ dE cm r(E cm ,T) , 

J Z7T Z TUj/^ J 

where r(E cm ,T) is the (unnormalized) collision spectrum for a given temperature T. Con- 
sequently, the average collision energy is given by 

< E cm {T) >= — J dE cm E cm r(E cm ,T) . 



3 It is conceivable that a nonvanishing parton content of the yet-to-be-formed hadrons leads to a 
finite dissociation cross section before Tp. Dimuon production has been shown to be enhanced by 
partonic interactions of preformed mesons in a similar context [18|. 
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III. RESULTS 



We first study very central Pb+Pb reactions (6 = 0) at _E lab = 160 GeV with the 
microscopic nonequilibrium model. Figure [I] shows the J / ip dissociation rates as functions of 
time. Comover dissociation processes occur most frequently about 1 fm/c after the nuclear 
reaction begins. At that time, the total dissociation rate reaches a value of more than 
0.1 c/fm. However, it is dropping rapidly, an order of magnitude within the first 10 fm/c 
of the reaction. The individual J/ipir and J/tp p interaction rates are plotted along with 
the total comover absorption rates. It is clear that 7r's and p's alone are not responsible 
for J /ip dissociation by comovers due to the large number of transient meson renonances 
present in the microscopic model. The string fragmentation scheme used in UrQMD |15| 
determines the relative abundancies of these states. This choice is far from unique since 
the model can be compared to data only for the hadrons observable in the final state. In 
UrQMD, all states of eight different meson multiplets ( J PC = 0"+, 0++, 1"", 1++, 2++ 
plus two excited 1 multiplets) may be populated in a string fragmentation process. The 
probability to form a meson from one of these multiplets is chosen to be proportional to the 
spin degeneracy and inversely proportional to the average mass of the multiplet. 

Table | shows the relative importance of the mesons which contribute most to J/ip dis- 
sociation. Scatterings of the J/^'s with 7r's and p's are indeed the dominant dissociation 
processes. However, together they are responsible for only 37% of the total comover absorp- 
tion. Twenty channels are left out of Table |, each of them contributing about 1% or less. 
However, they account for more than 15% of the total absorption. 

It is important to remember that the same dissociation cross section is used for all J/ip- 
meson interactions, most probably a crude approximation. For a better estimate of the J/ ip 
dissociation rates, one would need to calculate the cross sections in all possible channels. 
Of course, even if this whole set of cross sections was completely determined it is uncertain 
that the hadronic transport model describes all the details of rescattering correctly. In any 
case, it is questionable if a model of J/^-comover absorption which includes only the light 
mesons can be used for reasonable quantitative predictions. 

Figure ^ shows the average J/^-meson collision energies in central Pb(160 GeV)+Pb 
reactions as functions of time for J/ipn and J/ip p interactions. As one expects, the average 
collision energies are highest in the early stage of the reaction where the rates are also at 
their maximum. 

We now compare the microscopic nonequilibrium scenario with a thermal meson gas. 
Figure shows the rates and the averaged thermal dissociation energies in an ideal gas 
of 7r's and p's as functions of temperature. The rates in the later stage of the reaction, 
shown in Fig. [TJ, when t > 2 fm/c, correspond to a temperature of T = 140 ± 20 MeV in 
the equilibrated system. The temperatures one would deduce from the collision energies of 
Fig. |2] roughly agree with these values. Thus, the concept of a thermal hadron gas may 
be approximately valid in the later stage of a nuclear reaction. The survival probability of 
charmonium states, however, is to a large extent determined by the early reaction dynamics, 
t ~ 1 fm/c. Here, the rates roughly correspond to a thermal 7r and p gas with T « 220 MeV. 
In contrast, the average collision energies at this time are much higher than one could ever 
expect in a thermal hadron gas. As can be read off from Fig. ^| and Fig. |3], the collision 
energies would correspond to temperatures greater than 800 MeV. We can therefore conclude 
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that a major part of the J /if) dissociation processes cannot reasonably be approximated by 
a thermal scenario. Much higher collision energies are observed than can be expected in a 
thermal system indicating that the underlying tt and p momentum distributions are out of 
equilibrium and the probability of high relative velocities is enhanced. 

The nonequilibrium character of J/ip rescattering is also reflected by the time integrated 
collision energy spectrum. The prediction of the nonequilibrium model is shown in Fig. |], 
together with the collision spectra of a thermal gas at T = 200 and 300 MeV. Obviously, 
in the microscopic simulation collisions with high center of mass energy are much more 
probable than in a thermal hadron gas at reasonable temperatures. 

As mentioned above, we have used fixed dissociation cross sections, (Tj/^m = 2.41 mb, in 
this study. However, absorption by light mesons has recently been calculated in the frame- 
work of a meson exchange model p] where the cross sections are strongly energy dependent. 
Thus, at a fixed comover density, the assumption of thermal velocity distributions would 
considerably underestimate the absorption rates. At E cm = 5 GeV, the energy dependent 
cross sections of Ref. || are aj/^ n = 2.3 mb and <Jj/^ p = 1.6 mb. As a result, using these 
cross sections || in the nonequilibrium UrQMD simulation leads to slightly lower maximum 
absorption rates at t ~ 1 fm/c than do the fixed cross sections. In the later stages of the 
reaction, the energy dependent absorption rates are more strongly suppressed since most of 
the collisions occur closer to threshold, as shown in Fig. |2|. Typically, for t > 3 fm/c, the 
collision energies are E cm « 4 GeV which implies oj/^ ~ crj/^ p ~ 0.3 mb, considerably 
smaller than the cross sections when t ~ 1 fm/c. 

Certainly, different parametrizations of the energy and eigentime dependencies for all 
possible comover dissociation cross sections would yield different integrated survival probil- 
ities. In this work we explore only the influence of the comover dynamics on the absorption 
rates. Comparison to experimental data with a refined prescription of charmonium produc- 
tion and absorption in the microscopic framework will be presented elsewhere |14| . 



IV. CONCLUSION 

We have presented a microscopic model of J/ ip production and absorption in relativistic 
heavy ion collisions. Averaged dissociation rates and collision energies have been analyzed 
as functions of time. In this model, scattering of J/^'s with 7r's and p's account for only 
about 37% of the total absorption by comovers. The J/ip absorption by comovers is strongly 
time dependent with dissociation rates peaking at t ~ 1 fm/c. The collision energies at that 
time are much higher than can be expected in a thermal scenario. 

An equilibrated gas of light mesons seems to be too simple an approximation to the actual 
hadronic environment in which J/^'s interact. In order to reliably calculate the survival 
probabilities of charmonium states in heavy ion collisions, one must employ more realistic 
simulations of the meson dynamics. Many different reaction channels contribute to the total 
absorption. The corresponding dissociation cross sections and their energy dependencies 
need to be explored further with the help of phenomenological models. 
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TABLES 



TABLE I. Contributions of different meson species to the total J/ip comover absorption ac- 



cording to the UrQMD simulation. Only the most dominant channels are shown. 



Species 


Percentage of comover absorption 


7T 


18 5 ± 3 


n 
r 


18.4 ± 0.3 


K 


10.7 ±0.3 


K*(892) 


7.4 ±0.2 


V 


6.4 ±0.2 




6.4 ±0.2 


a 2 (1320) 


4.8 ±0.2 


ai(1260) 


4.3 ±0.2 


61 (1235) 


3.8 ±0.2 


a (980) 


3.0 ±0.1 


Sum: 


83.7 ±0.7 
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FIGURES 



all 




t (fm/c) 

FIG. 1. J/ip dissociation rates in central Pb(160 GeV)+Pb reactions according to the UrQMD 
calculation as functions of time. Shown are the rates for all J /ip- meson interactions (full line), 
those for J/ipir (dotted line) and those for J/ip p (dashed line). A fixed dissociation cross section 
of ojj^M = 2-41 mb is used for all mesons. 
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FIG. 2. Average J/^-meson collision energies in central Pb(160 GeV)+Pb reactions as a func- 
tion of time for J/tp-K (dotted line) and J/ip p interactions (dashed line) . A fixed dissociation cross 
section of cfj/^m = 2-41 mb is used. 
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FIG. 3. Left: J/ip dissociation rates in an ideal gas of pions and p's as a function of temperature. 
A fixed dissociation cross section of cfj/^m = 2-41 mb is used. Right: Average collision energies in 
the same system as a function of temperature. 
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FIG. 4. Left: J ftp collision spectrum in central Pb(160 GeV)+Pb reactions according to the 
UrQMD calculation (normalized to one collision). Shown are the energy spectra for J/tpn inter- 
actions (dotted line) and for J/ip p interactions (dashed line) . A fixed dissociation cross section of 
&j/ipM = 2.41 mb is used. Right: J/ijj collision spectra in an ideal gas of 7r's and p's for T = 200 
and 300 MeV, normalized to one collision. A fixed dissociation cross section of ojj^m = 2.41 mb 
is used. Note that the scales are identical in both plots. 
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